Roughness exponents and grain shapes 
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In surfaces with grainy features, the local roughness w shows a crossover at a characteristic length 
r c , with roughness exponent changing from ai ~ 1 to a smaller Q2.. The grain shape, the choice of 
w or height-height correlation function (HHCF) C, and the procedure to calculate root mean-square 
averages are shown to have remarkable effects on a.%. With grains of pyramidal shape, ai can be as 
low as 0.71, which is much lower than the previous prediction 0.85 for rounded grains. The same 
crossover is observed in the HHCF, but with initial exponent xi ~ 0.5 for flat grains, while for some 
conical grains it may increase to \\ w 0.7. The universality class of the growth process determines 
the exponents 012 = X2 after the crossover, but has no effect on the initial exponents qi and xi> 
supporting the geometric interpretation of their values. For all grain shapes and different definitions 
of surface roughness or HHCF, we still observe that the crossover length r c is an accurate estimate 
of the grain size. The exponents obtained in several recent experimental works on different materials 
are explained by those models, with some surface images qualitatively similar to our model films. 

PACS numbers: 68.35.Ct, 68.55.Jk, 81.15.Aa, 05.40.-a 
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I. INTRODUCTION 



Scaling properties of the local surface roughness w and 
of the height-height correlation function (HHCF) C are 
very useful to understand the growth dynamics of thin 
films and other deposits [IH3- The usual approach is to 
measure exponents from plots of w or C as a function 
of the box size r (roughness exponent) or time t (growth 
exponent) and to relate their values to some universal- 
ity class of growth [l|. However, a very small number 
of systems exibit simple scaling features to match those 
theories. For instance, the presence of grains in the film 
surface leads to a crossover between two regimes where w 
increases with r with different roughness exponents, ai 
and Q!2, as illustrated in Fig. 1 |3-ll0j. For the HHCF, the 
same crossover occurs with exponents xi an d X2-- Simi- 
lar crossover is observed in other systems, such as fresh 
snow on the ground and pyroclastic deposits on volcanic 
surfaces (lHll^. 

In Ref.pjl], the crossover with ai « 1 was shown to 
be a geometric effect of the grainy surface structure and 
of the gliding box method (analogous result is obtained 
with the box counting method). It was also shown that 
the crossover took place when r was close to the average 
grain size. If the grain surface is flat, a.% is very close to 
1, while for rounded grains it decreases to values close to 
0.85 [Hj]. These results match those of a large number of 
experimental works [3-iloj. However, other experimental 
works show film surfaces with grainy structure, the same 
crossover in roughness sca ling or HHCF, but with much 
smaller exponents a\ |14l - |22j . The usual interpretation 



for those exponents is that small scale surface features are 
determined by a different growth dynamics. Indeed, even 
the crossover with a\ « 1 was already interpreted as an 
anomalous scaling, with a± being called local roughness 
exponent (denoted a/ oc ) and a-i called global roughness 
exponent. For these reasons, in many systems it is still 
unclear whether a crossover similar to that in Fig. 1 
should be interpreted as a purely geometric effect or as 
a consequence of a competitive growth dynamics. 

Here we study several growth models with grainy sur- 
face features to show the possible effects of the grain 
shape, of the method of calculation of averages of squared 
quantities, of the working quantity (w or C) and of the 
universality class of the growth process. For all growth 
models, grain shapes, and methods of analysis, we ob- 
serve crossovers at box sizes very close to the average 
grain size. We also show that a very broad range of a.\ 
can be found, depending on the grain shape and the work- 
ing quantity, but independently of the universality class 
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FIG. 1: Typical behavior of the local roughness as a function 
of box size in grainy surfaces. 
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of growth, which determines only the value of ui. Sim- 
ilar conclusions are obtained for the exponents xi an d 
Xi- The comparison with experimental works with sev- 
eral materials and deposition methods gives additional 
support to the geometric interpretation of the crossover 
in those systems. 

The rest of this work is organized as follows. In Sec. 
im we define average quantities and present the growth 
models. In Sec. IIII1 we recall the results of some exactly 
solvable models with grains at the surface, which explain 
the crossover with a% « 1 and Xi ~ 0.5 (with the usual 
definition of the HHCF). In Sec. HVI we analyze the effects 
of the grain shape, particularly some very sharp grains, 
considering models in different universality classes. In 
Sec. El we show the applications of our approach to real 
films. In Sec. I VII we summarize our results and present 
our conclusions. 



II. DEFINITION OF AVERAGE QUANTITIES 
AND MODELS 

First we define the average quantities analyzed in this 
work. 

The surface roughness in square boxes of size r at time 
t is usually defined as 

w(r,t)^((h~h/ /2 y (I) 

The overbars in Eq. ([T]) denote averages of the height 
h inside a given box position (spatial average) and the 
angular brackets represent the configurational average as 
the box scans the whole surface of a deposit. This is 
called gliding box method, in which the scanning box 
moves one pixel each time it performs a new spatial av- 
erage. In box counting methods, the surface is divided 
in nonintersecting boxes for the configurational average. 
Alternatively, some authors define the roughness as 



w '=/(h-V) 2 \ , (2) 

i. e. they calculate the configurational average of the 
square height fluctuation and take the square root of that 
average. 

When several images of a deposit are available, or sev- 
eral configurations are grown with the same model, these 
different samples also contribute to the above configura- 
tional averages. 

For window sizes below the grain size, the roughness 
scales as 

w{r,t)~r a \ (3) 

which defines the initial roughness exponent u\ (Fig. 1). 

The height-height correlation function (HHCF) at dis- 
tance r and time t is usually defined as 

C(r,t) = ([h(r +r,t)-h(r ,t)} 2 ) 1/2 , (4) 



with configurational averages taken over all different ini- 
tial positions ro. Alternatively, it can be defined as 

C'(r,t) = (\h(r + r,t)-h(r ,t)\), (5) 

which corresponds to an interchange of the configura- 
tional average and the calculation of the square root in 
Eq. (UJ. In this sense, the calculation of C{r,t) parallels 
that of it/, while the calculation of C'(r,t) parallels that 
of w. 

For window sizes below the grain size, the HHCF scales 

as 

C(r,t)~r Xl , (6) 

which defines the initial roughness exponent xi f° r that 
function. 

For window sizes much larger than the grain size (i. e. 
r r c - see Fig. 1), a surface obeying normal scaling has 
w ~ r Q2 and C ~ r X2 , with «2 = Xi- The quantities w' 
and C obey the same scaling. Those exponents are rep- 
resentative of the large lengthscale kinetics governing the 
growth process. Typical examples of growth kinetics are 
those of Edwards- Wilkinson (EW) [23[ , of Kardar-Parisi- 
Zhang (KPZ) [24|, and the diffusion-dominated ones, lin- 
ear (Mullins- Herring - M H) |25I| or nonlinear (Villain-Lai- 
Das Sarma - VLDS) ^,[23^ 

Now we present the models for growth of thin films 
with grains at the surface. 

Intrinsic corrections to scaling for large r and large t 
should be avoided in those models, so that any crossover 
is solely due to the grainy structure. This request ex- 
cludes the grain deposition models introduced in Ref. 
[13} and related ballistic- like models [H, [2^] because they 
have remarkable scaling corrections. On the other hand, 
some models with smooth surfaces and particle enlarge- 
ment presented in fiefs. @, HH satisfy that condition. 
They are described below. 

The first model has KPZ kinetics. The first step is to 
grow a deposit with cubic particles of unit size following 
the rules of the restricted solid-on-solid (RSOS) model: 
the aggregation of the incident particle is accepted only 
if the height differences of nearest neighbors are always 
or 1 (otherwise the aggregation attempt is rejected) 
[30] . We recall that a 2 = Xi ~ 0.39 for the KPZ class in 
two-dimensional substrates |3lj |. 

The second model has VLDS kinetics. The initial de- 
posit is grown with the rules of the conserved RSOS 
model, where the incident particle executes a random 
walk between neighboring columns until finding a col- 
umn where it can aggreg ate respecting the conditions on 
height differences [32|3J]. We recall that a-i = Xi ~ 0.67 
for the VLDS class in two-dimensional substrates [331 l34j | . 

After growing the initial deposit, with KPZ or VLDS 
model, the size of each particle is enlarged by a factor I, 
i. e. each particle is transformed in a cubic grain of side 
I. Most of our simulations are performed with Z = 32. 
The final step is replacing the top cub grains (surface 
grains) by rounded or sharp structures. Three shapes 
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are used: semi-ellipsoids of horizontal radius ly/2/2 and 
vertical radius h, cones with that radius and height h, and 
pyramids of square basis of side I and height h. They are 
illustrated in Fig. [5] Several values of h are considered 
for each shape, typically between I and 3/. 




FIG. 2: Shapes of surface grains after enlargement of the 
original deposits: a) semi-elliptical, b) conical, c) pyramidal. 
Semi-ellipsoids and cones are cut at the sides so that their 
basis are squares of side I that fit the shape of the cubic grain 
at their bottom. This is the reason for radius lV2/2 of their 
basis. 

In the scaling of w or C, the role of the height h is 
measured relatively to the height of the surface steps, 
which is I. The horizontal scaling factor is also I for the 
cubic grains, but this is not important for the scaling 
exponents. For instance, if the grains were constructed 
with the shape of paralelepids of height I and horizontal 
sides Zy, the scaling exponents would not change. Thus, 
the aspect ratio of the grains considered here is not a 
limitation of the model. 

The simulations of the KPZ and VLDS models were 
performed in square substrates (three-dimensional de- 
posits) of lateral size L = 128 at times of order 10 4 . For 
the RSOS model, it corresponds to approximately 5 x 10 3 
layers of unit size particles; for the CRSOS model, cor- 
responds to 10 4 layers. After replacement of the original 
particles by grains of size I = 32, the deposits have lat- 
eral size 4096. Simulations in smaller sizes (L = 64 and 
L = 32 for the original models) and different grain size 
(/ = 16) give similar results for all exponents, indicating 
that finite-size and finite-time effects are negligible. 



III. THEORETICAL PREDICTIONS FOR ai 

AND xi 

As the scanning box glides along the surface, it fre- 
quently encloses high surface steps created between 
neighboring grains. These are the box positions where 
the largest height fluctuations are encountered, thus they 
give the main contribution to the roughness (Eqs. [1] or 
[2]). If the box has size r (i. e. r pixels in each direction), 
then the number of box positions that involve each high 
step is proportional to r. Thus, the configurational aver- 
age of Eq. |T]) gives roughness w proportional to r. This 
gives o.\ — 1, as explained in Ref. [13j and confirmed by 
simulations of several models. 

When Eq. @ is used, w' is a configurational average. 
The main contribution to that average also comes from 
box positions enclosing high surface steps, thus, that av- 
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FIG. 3: Scaling with the window size r of data for the KPZ 
model with cubic grains: a) w (blue squares) and w' (red 
triangles); and b) C' (blue squares) and C (red triangles). 



erage is proportional to r. This gives w' proportional to 
r 1 / 2 , i. e., ati = 1/2. 

These results are confirmed by our simulations of the 
RSOS model with cubic grains, as shown in Fig. 3a. It 
clearly shows the remarkable difference in the scaling of 
w and w' for box sizes smaller than the grain size, while 
the same exponent «2 after the crossover represents the 
universality class of the process. 

Similar situation is observed with the HHCF. Again 
the main contribution for the configurational average 
comes from box positions which involve high surface 
steps, thus this average is proportional to r. With the 
most used definition of that function (Eq. @|, we have 
C(r, t) proportional to r 1 / 2 , thus the crossover takes place 
with xi = 1/2- Instead, if the scaling of C'(r,t) is ana- 
lyzed, we expect xi = 1- 

Simulations of the RSOS model with cubic grains show 
the predicted crossover, as illustrated in Fig. 3b. The 
exponent \i is very close to 1/2 for C(r,t) and slightly 
below 1 for C'(r,t). Again, the universal exponent \i is 
obtained after the crossover; as expected, oli ~ X2- 

With the usual definitions of surface roughness (w - 
Eq. [J) and HHCF (C - Eq. gj), the roughness expo- 
nents measured before the crossover (at, Xi) are differ- 
ent. This contrasts with the expected universality after 
the crossover (ct2 ~ X2)- Our analysis show that those 



4 



discrepancies are effects of the grainy morphology and 
the calculation method, in particular the order of calcu- 
lation of square root and configurational average in Eqs. 
© and Q. 



IV. EFFECTS OF GRAIN SHAPE 

Rounding of the surface grains may lead to a\ between 
0.85 and 1, as shown in Ref. [13j. However, the replace- 
ment of the cubic grains by the rounded or sharp struc- 
tures in Fig. 2, with h > I, leads to much more drastic 
changes in the initial exponents of w(r,t) and C(r,t). 

This result is illustrated in Figs. 4a and 4b for films 
grown with the KPZ model and pyramidal grains of 
height h = 64: ol\ decreases to 0.71 and \l increases 
to 0.61 (while ct2 ~ 0.39). Fig. 4b also shows the forma- 
tion of a plateau in the HHCF before the second scaling 
regime, which is characteristic of all sharp grains with 
large heights. 

In Table I, we show the values of u\ and xi obtained for 
semi-elliptical, conical and pyramidal grains with several 
heights. 

A remarkable result is that films grown with the VLDS 
model have the same exponents ot\, Xi U P to the sec- 
ond decimal place, despite the significant change in the 





FIG. 4: a) Roughness (w) and b) HHCF (C) as a function of 
the window size r, for the KPZ model with pyramidal grains 
of heigth h = 64. 



FIG. 5: a) Roughness (w) and b) HHCF (C) as a function of 
the window size r, for the VLDS model with conical grains of 
heigth h = 32. 



asymptotic roughness exponent {a.% ~ 0.67). In Figs. 
5a and 5b, we show results for conic grains with height 
h = 32, which give a x = 0.809 and xi = 0.539. Those 
values are close to the KPZ values shown in Table I for 
the same grains. The main differences from the models 
with KPZ scaling are that the change in the slope of the 
roughness plot is smaller and there is a slope increase in 
the HHCF plot when passing from the first to the second 
scaling regime. 

Table I shows that ai is much smaller than the limit 
0.85 obtained in previous work [l3j for many grain 
shapes, particularly for sharp conic and pyramidal grains. 
With the structures studied here, the lower limit is close 
to 0.71, obtained with pyramidal grains. For h < 31, the 
general trend is that the increase of h leads to decrease 
of a.\. For larger h (not shown in Table I), a very slow 
increase of a\ towards 1 is observed. 

The relative changes in xi are much larger, attaining 
almost 40% for conic grains with h = 31 (see Table I). 
Indeed, this is the grain shape that provides higher devi- 
ations from the flat grain value xi = 0.5. A monotonic 
increase of xi is observed when taller grains are studied. 

The above results show that sharp grain shapes bring 
closer the exponents aj and Xij m contrast with the very 
different values for flat grains (1 and 0.5, respectively - 
Sec. IIIip . In some cases, they are surprisingly close; 
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h 32 64 96 





0.826 


0.773 


0.763 


xf E 


0.504 


0.551 


0.589 


a? 


0.806 


0.768 


0.768 


X? 


0.539 


0.633 


0.694 


a? 


0.755 


0.710 


0.708 


p 
Xi 


0.535 


0.606 


0.645 



TABLE I: Exponents obtained from w (ai) and C (xi) in 
KPZ films with semi-elliptical (SE), conical (C) and piramidal 
(P) grains. 



h 32 64 96 





0.523 


0.576 


0.621 


x! E 


0.754 


0.711 


0.705 


a? 


0.554 


0.650 


0.719 


X? 


0.712 


0.693 


0.710 


a? 


0.556 


0.638 


0.687 


xf 


0.664 


0.633 


0.640 



TABLE II: Exponents obtained from to' (ai) and C' (xi) in 
KPZ films with semi-elliptical (SE), conical (C) and piramidal 
(P) grains. 



for instance, they differ only 10% for conic grains with 
h = 31. 

In Table II, we show exponents «i and Xi obtained 
from the scaling of w'(r,t) and C'(r,t). They should be 
compared with the respective flat grain values 0.5 and 1. 

Comparison of results in Tables I and II show that 
sharp grain shapes also bring closer the values of a\ mea- 
sured from w and w' scaling, which are very different for 
flat grains (1 and 0.5, respectively - Sec. Illlj) . It is par- 
ticularly interesting to observe that ai differs only 3% 
when calculated from w or w' in films with pyramidal 
shapes with h — 31. These values may be incorrectly in- 
terpreted as true roughness exponents because the same 
Q!2 is expected for w and w' . This type of erroneous in- 
terpretation can be avoided if one accounts for the effects 
of a wide range of grain shapes and sizes and investigates 
other quantities, such as HHCF. 

The crossover size r c is defined at the intersection of 



the linear fits of the initial regime and the second scaling 
regime of roughness or HHCF, as illustrated in Fig. 1.. 
Despite the wide range of values of a\ and \i shown in 
Tables I and II, a remarkable result is that r c is always 
very close to the grain size I, for KPZ and VLDS models. 
Using I = 32, our estimates range between r c — 30 and 
r c = 34, which corresponds to a maximum difference of 
7%. Consequently, r c can always be used as a reliable 
estimate of the grain size. 

In the above models, we considered surfaces with uni- 
form grain height. However, we also analyzed the effect 
of distributions of grain heights, since this is the situation 
in real surfaces. In all cases, we observe that the expo- 
nents a\ and xi are near the averages of those obtained 
with a single value of grain height. 

An example of a film surface with such random grain 
distribution is shown in Fig. 6: 1/4 of the grains are flat 
and 3/4 have pyramidal shape, with equally distributed 
heights h — 32, h — 64, and h = 96. For that surface, we 
obtain a\ = 0.742 and xi = 0.554, which is close to the 
average of the results in Table I for those shapes. 




FIG. 6: Film surface with 1/4 of the grains flat and 3/4 pyra- 
midal with heights h — 32, h — 64, and h = 96 equally 
distributed. 



V. COMPARISON WITH EXPERIMENTAL 
RESULTS 

In the experimental works discussed below, the expo- 
nent a\ defined here is frequently named local roughness 
exponent a/ oc , as a reference to the small lengthscale be- 
havior and/or to a possible anomalous scaling. 

Several experimental works have already shown the 
crossover of Fig. 1 with a ~ 1, which is explained by 
the growth models with fiat or slightly rounded grains 
[l3l |. Among those works, we highlight the study of rf 
sputtered LiCoO x films by Kleinke et al Q, which gives 
0.91 < ot\ < 0.95; the spray pyrolysis growth of ZnO 
films by Ebothe et al @, which gives 0.94 < a x < 0.97 
for high flow rates; the electrodeposition of cooper by 
Mendez et al § and of gold by Vazquez et al Q , which 
give ai = 0.87 ± 0.06 and a x = 0.90 ± 0.06, respec- 
tively; the electrochemical roughening of silver electrodes 
by Otsuka and Iwasaki Q, which gives ot\ between 0.95 
and 0.98; and the pulsed laser deposition of La modied- 
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PbTiO^ films of Vasco et al which have ol\ = 1. 

However, many works show the same crossover with 
exponents a% between 0.7 and 0.85, and surface images 
confirm the presence of grains of approximately conic or 
pyramidal shape, much higher than the steps between 
neighboring grains. These features are observed in films 
of various materials and substrates, deposited with dif- 
ferent techniques. This justifies our approach with geo- 
metrical models, independently of the particular growth 
dynamics. 

Among the applications to inorganic materials, we find 
some vapor deposited gold films by Vazquez et al, which 
have ax w 0.83 - see Fig. lc and Fig. 3 of Ref. QJ]. 
One of the niquel oxide film samples deposited by sput- 
tering in Ref. [15[ have a.\ — 0.70, and the AFM image 
show the qualitative features of our models with sharp 
grainy structure. Nearly the same exponent {a\ — 0.71) 
is obtained with Ni films electrodeposited on indium tin 
oxide substrates in Refs. 0,1131 • Several Ni ~ Zn alloy 
films of Ref. [l8[ show the crossover in roughness scaling, 
with most estimates of a\ in the range [0.80,0.83]. This 
is consistent with our models of semi-elliptical grains of 
lower h, and the images actually show a smooth grain 
morphology. 

The same features are also observed in organic ma- 
terials. Films formed with bilayers of poly(allylaminc 
hydrochloride) and a side-chain-substituted azobenzene 
copolymer (Ma-co-DR13), after deposition of 10 or 20 
bilayers, show grains with a broad size distribution, and 
the initial roughness exponents 0.81 and 0.79 [jjj]. AFM 
images of chemically deposited polyaniline thin films on 
glass substrates [2fJ] have similar features, but, as far 
as we know, roughness scaling was not studied with 
those images. The surface of LangmuirBlodgett films of 
polyaniline and a neutral biphosphinic ruthenium com- 
plex (Rupy) of Ref. [2l| also show those grainy features 
with some high peaks, and initial roughness exponents 
are in the range 0.66 < ot\ < 0.81 for thicknesses be- 
tween 1 and 21 layers. However, most estimates of ot\ 
are between 0.72 and 0.76 [2l|, in good agreement with 
our results for pyramidal grains. 

It is interesting to observe that some surface images 
shown in Refs. [HI, [Tt], fl9l - l2l| have features similar to 
the model illustration in Fig. 6 (in most cases without 
the flat grains). This comparison reinforces our interpre- 
tation of the exponents measured in those works. 

Similar results are obtained in etching of silicon sur- 
faces in Ref. [221 ]: ct\ is found between 0.70 and 0.87 
when the (111) surface is etched by an NaOH solution 
in contact with a non-saturated aqueous environment. 

It is also important to recall that there are systems 
with crossover in the roughness scaling which do not show 
the sharp grainy features of our models, and consequently 
deserve separate investigation. For instance, the images 
of another sample from Ref. 15| does not show those 



features, but the roughness shows a slow crossover with 
ai = 0.52. Pyroclastic deposits of Mt. Etna show rough- 
ness scaling crossover with a\ between 0.47 and 0.67, but 



the images do not support modeling by grainy structures 
[l2l].. There are also systems with sharp grainy struc- 
tures and small a\, such as some Ni films of Ref. 
(0.55 < a\ < 0.61), which also would deserve a separate 
investigation (those films have 0.12 < a2 < 0.22, which 
also cannot be easily explained with the well known ki- 
netic growth theories [l|). 

The crossover in HHCF scaling obtained in some sys- 
tems can also be related to our models. For instance, 
Manes et al [ll| used HHCF as a measure of fresh snow 
roughness and obtained xi between 0.58 and 0.62 in a 
set of five experiments. These values are consistent with 
our model with very high semi-ellipsoidal grains {h = 3/) 
or with conic or pyramidal grains with h = 21 or less. 

Again, there are also systems where a crossover of 
HHCF scaling is observed but whose images do not show 
the features of our models. An example is Ref. [35j |. 
where xi — 0-84 was obtained for paraphin films de- 
posited on stainless steel covered with amorphous car- 
bon. 

Results of the recent work on pentacene island growth 
on stepped oxide surfaces [36| can also be related to our 
models. First, for long lengths, the HHCF has exponent 
2x — 1, which is expected for height fluctuations domi- 
nated by the surface steps; indeed, arguments analogous 
to those for flat grains (Sec. IIIip give x = 1/2- How- 
ever, for small lengths, height fluctuations in the surface 
terraces (due to pentacene islands) lead to the increase 
of the HHCF exponent to the range [0.69,0.8]. Recent 
works showing evidence of anomalous scaling in organic 
and inorganic film surfaces also give estimates of HHCF 
exponents above 1 /2 at short lengthscales [37L [38| • Al- 
though both short and long range dynamics may be much 
more complex than in our models, a simple geometric in- 
terpretation of the short range exponents may also be 
considered due to the presence of grainy structures in 
the surface images. 



VI. CONCLUSION 

We extended the work on growth models with grainy 
surfaces to analyze the effects of the grain shape, of the 
method of calculation of averages of squared quantities, 
of the working quantity (roughness or HHCF) and of the 
universality class of the growth process. For all mod- 
els, grain shapes, and methods of analysis, we observe 
crossovers at box sizes very close to the average grain 
size. We also show that a very broad range of the initial 
exponent a\ is found for the roughness scaling, decreas- 
ing from 1 for flat grains to 0.71 for some sharp pyra- 
midal grains. The initial exponent xi of HHCF scaling 
increases from approximately 0.5 for flat grains to val- 
ues larger than 0.7 for sharp conic grains. Simulations 
of KPZ and VLDS models show that the universality 
class has no significant effect on the estimates of a\ and 
Xi- The range of a\ presented here explains results of 
some recent experimental works with different materials 
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and deposition methods. This gives additional support 
to the geometric interpretation of the crossover in rough- 
ness scaling for a variety of systems. 
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